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We address the problem of quantum process tomography with the preparators producing states 
correlated with the environmental degrees of freedom that play role in the system-environment 
interactions. We discuss the physical situations, in which the dynamics is described by nonlinear, or 
noncompletely positive transformations. In particular, we show that arbitrary mapping Qi n — > g ou t 
can be realized by using appropriate set of preparators and applying the unitary operation SWAP. 
The experimental "realization" of perfect NOT operation is presented. We address the problem 
of the verification of the compatibility of the preparator devices with the estimating process. The 
evolution map describing the dynamics in arbitrary time interval is known not to be completely 
positive, but still linear. The tomography and general properties of these maps are discussed. 



I. MOTIVATION 

The postulates of quantum theory require that the 
dynamics of isolated quantum systems is driven by 
Schrodinger equation [1, 2], i.e. for each time interval 
the evolution is described by a unitary transformation. 
However, for open quantum systems the situation is dif- 
ferent [3, 4] and under certain assumptions the evolution 
is described as an one-parametric sequence of completely 
positive tracepreserving linear maps (quantum channels) 
£ t - These mappings describe the state dynamics for ar- 
bitrary time interval (0, t), however for general time in- 
tervals (tijt?) the state transformations £t t ,t 2 '■ Qt t — * Qt 2 
do not necessarily possess the above property of complete 
positivity. The aim of this paper is to analyze the cases, 
in which the description of quantum dynamics is not com- 
pletely positive, or even not linear. One of the discussed 
problems will be the question of properties of the evolu- 
tion map £ti,ta fo r intermediate time intervals for general 
dynamics of open system governed by sequence £t- An 
important exception is if the sequence £t fulfills the semi- 
group property, i.e. £ t + s — £t£s for all t,s>0. In this 
case for each intermediate time interval the dynamics is 
linear and completely positive. 

The lack of complete positivity and linearity for state 
transformations is usually interpreted as an unphysical 
property, i.e. these operations cannot be physically re- 
alized. But still an optimal physical approximation of 
several physically impossible processes is of great impor- 
tance. Typical examples are quantum NOT operation 
[5], quantum copy operation [6, 7], etc. These processes 
violate the rules of quantum dynamics. However, we will 
see under which circumstances and in which sense even 
such unphysical transformations can be observed in our 
labs. 

Quantum process tomography is a particular goal of 
quantum experiments. Several strategies how to gather 
valuable experimental data for this task and methods 
how to correctly proceed such data are known [8-14]. 
The failure of the direct (inverse) estimation methods 
that could result in an unphysical map [9, 13], is usu- 



ally corrected by usage of more sophisticated statistical 
tools such as maximum likelihood [10, 11], or Bayesian 
statistics [12]. These methods are "forced" to lead to 
a correct quantum channel. We used to say that the 
failure of the direct estimation schemes follows from the 
finitencss of the measured statistical sample, i.e. the ob- 
served frequencies do not correspond to theoretically al- 
lowed probabilities and consequently, they do not corre- 
spond to some completely positive tracepreserving linear 
map. Here we shall address the question, under which 
circumstances such "unphysicality" can be due to imper- 
fect (potentially correlated) preparators. 



II. NONIDEAL PREPARATIONS 

The usual picture of open quantum system dynamics 
is based on three assumptions: i) the physical object un- 
der consideration is a part of some larger system that is 
isolated, i.e. its evolution is unitary, ii) initial state of 
the object and the environment is factorized, and iii) the 
state of the environment is independent of the state of the 
system. Under such conditions the resulting dynamics is 
completely positive and linear. The question is whether 
the unphysical maps obtained as a result of direct esti- 
mation can be interpreted in this picture provided that 
we relax the last two conditions, i.e. the initial state is 
potentially correlated, or the state of the environment 
depends on the system state, or both. This question, 
namely how the initial correlations affect the dynamics of 
open system, has been already studied by several authors 
[15-20]. Authors in [17, 19, 20] analyze this problem and 
propose the mathematical tools how to mathematically 
describe such extended evolution maps. 

A general state of bipartite system can be written in 
the following form 



Qab = Qa® Qb 
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efficients F jk = (Af®A%) eAB - (A/) eA (Af ) ea form the 
so-called correlation matrix. The evolution £ of the sub- 
system A is a composition of the following three maps: 
1) preparation map [21] V : Sa — » Sab (Sx stands 
for the set of quantum states of the system X) satis- 
fying the property TybVIqa] = QA, 2) isolated dynamics 
U : Sab — > Sab, i- e - = ^[2ab] = Uq A bU\ and 3) 
partial trace Tg : Sab — * Sa- The last two mappings are 
linear and completely positive. Moreover, both of these 
two features are preserved under the composition of map- 
pings. That is, the only source of "nonphysicality" is the 
preparation map V . One can show [15] that the linearity 
of the resulting dynamical map £ = Tb ° U o V requires 
V be of the form V \qa\ = QA <8> £b with £b arbitrary, but 
fixed. In such case the linearity and complete positivity 
of £ holds. 

In [20] authors studied different types of preparation 
maps and define the notion of an accessible map. The 
transformation is called accessible if it can be written 
as a composition of the preparation, unitary transforma- 
tion and partial trace [22] . If one allows arbitrary initial 
correlations, then the transformation is composed of two 
terms [16] 

Q ' A = ]T A^qA\ v + J2 1> 5>l A/ ® A£| M ) . (2.2) 

iiv jk u 

where the operators A^ v = {n\^/p^U\v) depends on qb, 
because are eigenvectors of the operator qb- That is, 
even if we put T = 0, the transformation g A — > q'a is still 
not necessarily described by some proper quantum chan- 
nel, because the choice of qb specifying the preparation 
map V can depend on g A - Consider an arbitrary state 
transformation Qi n — > g ou t- Let us define a preparation 
in the following way V [Qin] = Qin ® Qout- Next apply 
the SWAP operation (this is a unitary transformation) 
to obtain U S wAp(Qin ® Qout)U~l WAP = g out ® Q in . After 
performing the partial trace we obtain the required state 
transformation g in — > Q out . Moreover, we did not use any 
correlation (quantum, or classical) in our preparation at 
all. Of course, this construction is a bit artificial, but 
nevertheless it shows that arbitrary state transformation 
Qin — * Qout is in principle accessible, i.e. can be written as 
TboUoV. In order to avoid such "artificial" realizations 
of any map we need to pose some well-motivated physi- 
cal conditions. In the Ref.[17, 19, 20] the authors restrict 
themselves to linear preparation maps. In what follows 
we will analyze two experimental situations in which a 
preparation map is naturally defined and the extended 
dynamics can be studied. 

Before we get further let us mention one very impor- 
tant implication of the fact that arbitrary channel is 
accessible. In a sense this statement is very positive, 
because whenever we find out in our experiments some 
strong evidence that the dynamics is not linear, or not 
completely positive, we cannot automatically conclude 
that the quantum theory is not correct. The observed 
"unphysicality" can be still interpreted as the problem of 
devices called preparators that produce states correlated 



to degrees of freedom relevant for the subsequent system 
evolution. Without considering the dynamics the poten- 
tial correlations are not interesting and from kinematic 
point of view they are irrelevant. However, these dynami- 
cal aspects can be used to differentiate between otherwise 
kinematically equivalent preparators. From this point of 
view the "nonphysicality" means that the preparators are 
not independent of the environmental degrees of freedom 
that do take a role in the dynamics. 

As an example consider now the realization of the per- 
fect NOT operation realized on pure states \ip) — > \tp±)- 
Let us assume that the preparation of the spin-^ pure 
state is performed by a postselection after Stern-Gerlach 
measurement. In this way we can prepare any pure quan- 
tum state and mixtures can be obtained by mixing the 
pure state preparations. Kinematically this is completely 
correct preparation procedure of arbitrary qubit state. 
Imagine a situation that the spin is entangled with an- 
other spin such that together they are described by the 
singlet, i.e. |*) = ^(|V>)IV>-l) ~ \ip±)\ip))- Reading the 
outcome tells us perfectly which state \ip) we prepared, 
but the measurement affects also the state of the second 
spin, which is described by the state \ip±}- If the un- 
known device internally just swaps these two spins, than 
we find out that the device performs the transformation 
\tp) — > \tp±), i-e. the perfect NOT operation. Experi- 
menter using such preparations is not aware of the initial 
correlations and therefore he would conclude that the un- 
known device performs perfect quantum NOT operation. 
This conclusion is not wrong and indeed experimenter 
can prove that this device performs NOT operation, but 
only for specific preparation procedures. If he would use 
different state preparators (i.e. kinematically equivalent 
to the previous ones), he will very soon find some con- 
tradiction. 



III. PROCESS TOMOGRAPHY 

The "nonphysicality" is not an exception in process 
tomography using the direct estimation schemes [24]. 
In such schemes we usually measure a collection of as- 
signments Qj — > Qj for linearly independent states {Qj}- 
Using the linearity of quantum channels these assign- 
ments provide us with sufficient information to complete 
the reconstruction task. In other words: each state q 
can be written as a linear combination of states Qj, i.e. 
q = ^2j a j Qj ( a j ar e arbitrary) . Therefore the transfor- 
mation of q is determined by the set of measured assign- 
ments. However, quite often the resulting map is not 
completely positive. Even in cases when all Qj , g'j are 
proper quantum states. A reason could be really only 
the usage of finite data sample, i.e. the statistics is really 
small to conclude something about actual probabilities, 
mean values, and states. This line of arguments leads 
us to the usage of sophisticated statistical techniques 
(maximum likelihood, Bayesian approach, etc.). These 
statistical methods are modified for quantum tomogra- 
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phy purposes in a way that they are essentially forced to 
guarantee a physically valid result (a completely positive 
map in quantum process tomography) even in situations 
when the direct inverse schemes fail. However, as we have 
just explained, this lack of complete positivity can be due 
to presence of correlations in our preparators as well. If 
this is the case it is really not easy to say which of the 
preparators are not perfect, i.e. which of them represent 
the source of nonphysicality. It could happen that only 
one of them is imperfect, or all of them are imperfect. 

Without any insight into the physics behind the prepa- 
ration process, the measured data do not contain any in- 
formation about the origins and form of the preparation 
map. As we have seen there is always a trivial example 
using the SWAP operation that can be used to interpret 
arbitrary result. It is important to say here, that even if 
the "linearization" of assignments gives a correct quan- 
tum channel, it does not mean that our preparators are 
perfect. We should have in mind that the linearity is not 
tested, but only used as the theoretical tool to accom- 
plish the process estimation. In order to be sure about 
linearity one really needs to test the action of the chan- 
nel on preparators of any quantum state. In this sense 
the specification of the channel is always only a hypoth- 
esis and for specific (imperfect) preparators we can find 
the channel to be "unphysical" . Important point is that 
the usual notion of quantum channel has a good meaning 
only for properly prepared input states, i.e. for compat- 
ible preparators. 

We said that using only the measured data we have 
only very partial information about the real physics of 
the process. There are many possible unitary representa- 
tions of the observed assignments. One of the option is to 
say that a collection of used preparators was not perfect 
and conclude that the process estimation is not possible. 
Another approach is to apply some techniques of incom- 
plete process tomography [27] to estimate the unitary 
map acting on the system with Hilbert space Ti. ® H-env, 
where Tl env is arbitrary. However this is indeed a diffi- 
cult task, because we need to deal with data that do not 
contain complete information about the inputs as well as 
about the outputs of the channel, i.e. the assignments 
are not completely known. The method that can be used 
in such situations is called principle of maximum entropy 
[28]. However, these issues are beyond the scope of this 
paper. 



IV. COMPATIBILITY OF PREPARATOR AND 
TESTED PROCESS 

To be sure that the channel estimation gives a physical 
result one needs to use the collection of "good" prepara- 
tors producing linearly independent test states. "Good" 
in a sense that whatever degrees of freedom enters the 
preparation process, these are irrelevant for the channel 
realization. We will say that such preparators are com- 
patible with the channel realization. Let us note that 



the condition of producing a factorized state is not suffi- 
cient, i.e. even pure state preparators are not automati- 
cally free of imperfections. This follows from the example 
with SWAP operation, where no correlations in prepara- 
tor process are used at all, but any transformation can 
be realized. We see that important question is: how to 
test the quality of the preparator, or better to say, how to 
test the compatibility of the preparator and the quantum 
process? 

The motivation for the scheme we are going to use 
comes from the preparation process used in real exper- 
iments [25, 26, 29, 30]. In some cases the preparation 
of different states is done by exploiting quantum pro- 
cessing, i.e. transforming the known state by a known 
transformation to obtain a new state. In particular, let 
us assume we have a preparator that produces system in 
a state g. Applying unitary rotations Uj we are able to 
prepare states Qj — UjgUj that can be used to test the 
properties of an unknown quantum channel. Our aim is 
to test the compatibility of the original preparator and 
some unknown device (black box). Except the case of 
g = ^1 the unitary processing produce sufficiently many 
linearly independent states to perform the complete pro- 
cess tomography. This procedure, of course, requires per- 
fect realization and control of unitary transformations 
Uj. Moreover, these unitaries must be already "compat- 
ible" with the preparator. In some sense we are cheating 
a bit here, because we are going to test the compatibility 
of the preparator with a given device and we already as- 
sume we have devices (performing unitaries) compatible 
with the preparator. However, as we said this is quite 
usual procedure how to prepare different states in many 
experiments. Therefore, let us assume that we indeed 
have such compatible devices. Then the described setting 
can be used to test the quality of the single preparator 
with respect to the realization of the unknown channel £. 
Important point is that these "preparing operations" are 
independent of the original preparator, i.e. they do not 
introduce the "unphysicality" and the only source of "un- 
physicality" is the preparator of g. Let us note that there 
is no need for preparing operations to be unitary, but in 
general, they must be linear and completely positive in 
order not to introduce extra sources of "unphysicality" . 

In what follows we will study the properties of the 
preparation map given the described model of the ex- 
periment consisting of a preparator, an unknown black 
box £ and a set of preparing operations {$} (not neces- 
sarily unitary ones). Consider that the original prepara- 
tor (the one we want to test with respect to the chan- 
nel action) produces states ui = g ® gs + L, where 
L = Y^jk Fjk-hj ® Afe stands for the correlations. After 
fixing the set of preparing operations {$} the preparation 
V is, in general, defined as follows 

V : £>$ ^ = <f>®l[Lu} = <g> g B +T' (4.1) 

with r = $ ®j[r] = £ . fc r^A, ® A fe (V jk = ^ jl v lk , 

= Tr(Aj$[A/s])) and g$ = $[p]. The possible prepa- 
ration mappings V are represented by subsets of com- 
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pletely positive tracepreserving linear maps <I> that trans- 
form the given state g into arbitrary state in the domain 
of V . The simple example with the SWAP gate is ex- 
cluded/trivial in this case, because it generates only con- 
tractions into the fixed state qb- Moreover, the correla- 
tions induced in the preparation process are in some sense 
fixed by the state ui, i.e. by the original state preparator 
that we are testing (together with the channel reconstruc- 
tion). 

The unitary evolution U = v A^ v ® induces 
the state transformation 



0$ — > 
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where we used that qb = J2» ^A v ){ v \i i- e - {W)} 
are eigenvectors of qb- Since qb is fixed for any op- 
eration <!> we see that the first term is independent 
of <i>, i.e. it is independent of the preparation map 
V. Consequently, the unknown transformation £ de- 
scribing the device can be written as the sum of lin- 
ear completely positive map T {T\g\ = F^gF^ 
with F^ v — y/p^A^) and some traceless operator [20] 

s e = Ei,i, ft ,/ r ;-J A tU'V A i^. Le - Q - Q' = 
£[g] = F[g] + S e . First part is irrelevant of the correla- 
tions, but the second part can be even nonlinear. The 
properties of S e depends completely on the choice of the 
set of preparing operations i.e. choice of preparation 
mapping V . In special cases it can be linear and not 
factorized, but then it cannot be defined on the whole 
state space. The linear case was analyzed in Ref.[20]), in 
which the sufficient conditions for its existence (in terms 
of properties of S e ) were formulated. The linearity of 
V corresponds to a specific choice of the set of prepar- 
ing operations {$} for a given ui, but the specification of 
particular conditions remain to be an open problem. 

Another open question is the characterization of all 
transformations (not only with linear preparation map 
V) that can be understood within this model. This is 
indeed a very interesting, but also very difficult prob- 
lem and we are not going to discuss this here. Our aim 
is to describe the idea how to test the compatibility of 
the original preparator and the action of an unknown 
device. The compatibility means that either the initial 
correlations vanish (r = 0), or the unitary transforma- 
tion generating the process dynamics is from the set of 
transformations U <E {Ua <8> Ub, Ua <8> ^b^swap}- The 
pure state preparators are specific examples of prepara- 
tors with vanishing initial correlations (r = 0), and using 
the described procedure the pure state preparators are 
always compatible with an arbitrary quantum process. 

For instance, consider the preparator of a single pure 
state described in the example of the perfect NOT real- 
ization. That is, consider initially the spin is maximally 
entangled with another spin and by measuring along the 
z direction we are preparing the states | |)- Instead of 
using different measurement apparatuses to generate ar- 
bitrary pure states, let us perform single qubit rotations 



of the state | |) to create arbitrary pure state As be- 
fore, the device just swaps the two spins, i.e. arbitrary 
is replaced by | {). Consequently, the transformation we 
obtain in the single point contraction of the Bloch sphere, 
i.e. linear and completely positive map £[g] = | j)U I- 
That is, the same black box can be described by differ- 
ent quantum processes depending on the properties of 
the preparation procedures. 

To test the dynamical compatibility of an unknown 
preparation device we suggest to exploit calibrated de- 
vices performing some well known quantum operations. 
In fact, this is nothing new, because the same procedure 
[Afc]„,y' A^AjAj^iis^scd in most of the experiments. Important point is 
that using such method the properties of the resulting 
preparation map V (and consequently £) strongly de- 
pends on the properties of the original preparator. If one 
observes some "unphysicality" of £ in such experimental 
setting then before applying "statistical corrections" one 
should verify the compatibility of the preparator. Per- 
forming the process tomography experiment the set of 
preparing operations should be chosen in a way that the 
generated set of states is sufficient for the process tomog- 
raphy, i.e. this set is finite. The goal of the verification is 
to find whether the correlation matrix Y vanishes, or not. 
There are, in principle, two strategies that can be com- 
bined. We can use different preparing operations $i,<E>2 
generating the same state g^ 1 = p$ 2 and see whether the 
channel action generates the same state. Alternatively, 
we can use additional preparing operations to verify the 
linearity. In order to see, whether the unphysicality is 
due to imperfections in the preparation process, a char- 
acterization of all preparation maps V in the described 
settings is needed. This characterization is an open prob- 
lem the solution of which is necessary if we want to be 
able to propose some universal verification and estima- 
tion strategies. 



V. EVOLUTION MAP FOR ARBITRARY TIME 
INTERVAL 

Let us assume that the time evolution of the system is 
described by a set of completely positive maps £ t induced 
by some underlying unitary dynamics U t of the system 
and its environment. This corresponds to a situation in 
which initially (time t = 0) the system and the environ- 
ment are factorized, i.e. ojo — Q ® £■ Only under such 
assumption the maps £t can be completely positive for 
all t and we have £ t = TrU t g <£> £f7 t . In this section we 
turn back to the original question posed at the beginning: 
what are the properties of the time evolution maps £t ± ,t 2 
(t 2 > ti) describing the dynamics during arbitrary time 
interval [ti,^]? 

A direct calculation gives us that 



(5.1) 



These maps are linear, tracepreserving and hermiticity 
preserving, and thus defined on any operator (quantum 
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state), but except g e S = £t 1 [S(H)] C S(H) they can 
transform quantum states into negative operators. Fol- 
lowing the notation of Refs. [17-20] we define the positiv- 
ity domain T> pos (£) of a linear map £ as a subspace of 
states on which £ is positive, i.e. V pos (£) ={g£ S(H) : 
£[q] > 0}. In fact, physically, only the action on such 
subset is important, because this is what we can really 
test in our experiments. Let us note that £^ corresponds 
to the mathematical inverse operation (not physical) and 
does not necessarily always exist. This means that the 
question about the form of £t lt t 2 does not make sense if 
t\ is a point in which the maps £ tl is not invertible. 

The described model illustrates another physical situa- 
tion, in which the evolution map extended to whole state 
space is not completely positive. But in this case it is still 
linear, what makes its description much simpler. More- 
over, since the evolution maps £t 1: t 2 preserves the trace 
and hermiticity, they can be expressed as a difference of 
two completely positive maps [17, 23] 

q d 2 

£t ut2 [Q] = E A iQ A ] E A M , (5.2) 
j=l j=q+i 

where the operators Aj can be chosen so that they 
form an orthogonal operator basis (TrAjAk = for 
j 7^ k) and d is the Hilbert space dimension of the sys- 
tem. The tracepreservation is reflected by the identity 

Sj=i A ) A i - J2j= q +i A ] A i = !■ 

The physics behind such form of noncomplete positiv- 
ity is simple and just reflects the fact that at time t\ 
the system is correlated to relevant degrees of freedom 
that affects the forthcoming evolution. In each time in- 
stance t the global state of the system plus environment 
is described by some u) t = UtUJoU-t- The preparation 
map in time t is determined by the choice of w = Q ® £■ 
In particular, V t [gt] = w t = U t Lu U- t = U t Vo[go\U- t = 
UtVofi-^gtW-t, i.e. P f =W t o? o £~\ That is, the 
evolution for time interval (t, t + 5t) can be written as 
follows 

£st[gt] = TT B [U 5t V t [gt}ui} = £ t+St ° £ t _1 [f?t] . (5.3) 

Let us note that we can generalize the whole setting by 
allowing arbitrary preparation map Pq, but we want to 
preserve the physical picture with the factorized prepa- 
ration and therefore we will restrict ourselves to linear 
and factorized initial preparations only. 

Before we get further let us note that the problem of 
quantum process tomography for linear noncompletely 
positive maps was discussed in dissertation of Anil Shaji 
[19]. To our knowledge this was the first attempt to un- 
derstand the observed data in more general settings than 
just in the framework of completely positive maps. In 
particular, he analyzes the problem from the mathemat- 
ical point of view. Since the maps of this form are linear, 
the reconstruction schemes based on complete data are 
the same independently whether the complete positiv- 
ity constraint is applied, or not. The inverse estimation 



schemes use just the linearity of the quantum evolution 
and the observed data are collected so that the result is 
represented by some linear transformation uniquely. The 
only and crucial question is what type of linear transfor- 
mation it is. According to usual model of open system 
dynamics we expect to obtain completely positive trans- 
formations in our experiments, but in reality this is not 
always the case [29, 30] . Of course, the origin of this phe- 
nomena is questionable, but correlations of preparation 
map within the discussed model provide one possible op- 
tion. The estimation schemes and algorithms must be 
modified accordingly in cases, when indirect statistical 
methods (such as maximum likelihood, or Bayesian ap- 
proach) are employed, or our information is incomplete. 
Those who are interested in the details of complete quan- 
tum process tomography for noncompletely positive, but 
linear maps we refer to [19]. In this paper we are propos- 
ing the corresponding physical situation specifying the 
conditions under which the linear noncompletely positive 
maps can be observed experimentally. 

Our interest is not only to perform the process tomog- 
raphy experiment and reconstruction, but also to under- 
stand (at least partially) the physics behind. As we said 
the evolutions derived from the unitary dynamics U t for 
a fixed time interval are linear, tracepreserving and also 
hermiticity preserving. Therefore, they are of the form 
as written in Eq.(5.2). We are interested in the inverse 
question, whether any such transformation £ can be un- 
derstand as subdynamics between two instants of time 
induced by some unitary dynamics. Or alternatively, un- 
der which circumstances the preparation mapping is lin- 
ear on subset of quantum states and whether these situ- 
ations can be always understand as part of the dynamics 
described by £ t derived from unitary dynamics U t . 

If we assume that the one-parametric family of uni- 
taries Ut is arbitrary (i.e. the generating Hamiltonian is 
highly time dependent) , then there are no constraints on 
the choice of the transformations £ t for different times. 
One can always define the generating unitary transfor- 
mations Ut so that Tb o Wj o p = £ t . Thus the ques- 
tion is, whether the following identity can be fulfilled 
£st = £t+8t £t 1 -> where on the left hand side we have 
arbitrary linear transformation given by Eq.(5.2) and on 
the right hand side we have two arbitrary completely pos- 
itive maps £\,£2- The inverse operation £^ 1 is linear, 
tracepreserving and hermiticity preserving as well, i.e. it 
serves as the potential source of noncomplete positivity. 
It is a well known fact that substraction of two com- 
pletely positive maps realizes arbitrary noncompletely 
positive linear map (Eq.(5.2)), but here the question is 
whether a similar property holds for the "division" of 
two completely positive maps, i.e. for the transforma- 
tion £ t+ st °£ t -1 . 

Consider now the following example. The transposi- 
tion £trans is probably the best known noncompletely 
positive linear map. Its action is defined as follows 
^translf 1 ] = g T and for qubit it is closely related to per- 
fect NOT operation. Let us assume that as the result of 
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the process tomography we obtain the NOT operation, 
£not[£>] = \ [px Q&X +c y Qv y +<^zQ (7 z~q)- Could it happen 
within the discussed framework? Both these maps are 
positive, i.e. the positivity domain equals to the whole 
state space. In our settings the positivity domain always 
corresponds to an image of the whole state space under 
some completely positive map £t- However, only for uni- 
tary transformations the image of the state space equals 
to its original, i.e. if £t is a linear completely positive 
map and £ t [S(H)] = S(H), then £ t is unitary. Therefore 
we have £ t = U. Because of the unitarity this process is 
invertible. Consequently, Et+st °£^ 1 = £t+st = Est 

is necessarily a completely positive map. But this is in 
contradiction with the fact that our reconstruction gives 
us a noncompletely positive linear map £not (or £trans)- 
It means that the NOT operation £not, or transposition 
£trans cannot be interpreted as an evolution map describ- 
ing the time dynamics between two instants of time gen- 
erated by a global unitary dynamics with initially fac- 
torized preparation map. In fact, the same conclusion 
holds for arbitrary positive (but not completely positive) 
linear map transforming pure states onto pure states, i.e. 
whenever the identity £[S(H)] — S(H.) holds. 

However, there is still an option how to "partially" per- 
form the perfect NOT operation in the given framework 
of intermediate dynamics. The depolarizing single qubit 
channels form a one-parametric family £{ x } '■ t ~ * xr, 
where r is the Bloch vector corresponding to a quantum 
state g = \{I + f ■ a). For x e [—1/3, 1] the transforma- 
tions £{ x } are completely positive and £not = £{ x =-i}- 
For the inverse operations we have = £{i/ x } an d 
for the composition £{ y ) ° £{ x } — £{x.y} for arbitrary 
real x,y. Using all these identities it simple to proof 
that £not = £{ x } ° £{- x } = £{x.(-i/x)} = ^{-i}- This 
identity makes sense only for —1/3 < x < 1/3, when 
both transformations are completely positive. Formally, 
the above calculation suggests that we are able to realize 
the perfect quantum NOT operation during the dynam- 
ics governed by one-parametric set of completely positive 
maps. But, the above decomposition possessed a physical 
meaning only for states from the subset £{ x y[S(H)], i.e. 
for states with Bloch vectors smaller than \x\ (|f| < \x\). 
The maximal set of states on which we are able to re- 
alize the perfect NOT operation (in the given model) is 
contained in the sphere with radius x = 1/3. The con- 
clusion is that the perfect NOT operation can be find 
out as a result of the process reconstruction. Moreover, 
it can be even understood as an intermediate dynamics, 
but it does not mean that the perfect NOT operation is 
indeed accomplished, because no intermediate dynamics 
can perform a perfect NOT operation for all set of states. 
Thus, process estimation of the quantum operation be- 
tween two instants of time could result in perfect NOT 
operation. But the perfect NOT operation is performed 
only on restricted set of states. 

Although the presented framework of open system dy- 
namics enables us to explain quite naturally the exper- 
imental evidence of linear noncompletely positive maps, 



the answer to the inverse question is open, i.e. whether 
all hermiticity preserving, tracepreserving and linear 
transformations can be interpreted within the described 
model, if we relax its physical validity for the whole set 
of states. The characterisation of those noncompletely 
positive maps that can be realized within the discussed 
physical model is an open question that indeed requires 
deeper investigation. 



VI. CONCLUSION 

In this paper we have analyzed the dynamics of open 
quantum systems beyond the complete positivity restric- 
tion and related consequences for the process tomogra- 
phy. The correlations can be detected if the direct esti- 
mation procedure gives physically invalid result and si- 
multaneously, all the experimental and statistical devia- 
tions can be excluded. The good news of our analysis is 
that any failure of the process estimation (i.e. "unphysi- 
cal" result) cannot be interpreted as the failure of quan- 
tum theory unless one can exclude the presence of initial 
correlations and dynamically incompatible preparators. 
However, the particular realization using the SWAP op- 
eration is quite artificial. As an example, we have de- 
scribed two different experimental realizations of the per- 
fect NOT gate, which is considered to be unphysical. The 
bad news is that although the initial correlations can be 
detected, the process tomography is very difficult and 
ambiguous, and the physical origin could be quite arti- 
ficial. However, the usage of statistical tools is justified 
only if we can safely exclude all such (artificial) possibil- 
ities, i.e. we have some restrictions and models on the 
form of possible preparator maps. 

We have discussed two very natural physical situations 
that can result in observation of initial correlations ef- 
fect. First of them is motivated by current experiments, 
in which experimenters typically use single state prepara- 
tor and other states are generated with the help of fur- 
ther processing, for instance by applying different unitary 
operations. In the second case the "unphysicality" is re- 
lated to the fact that the evolution map describing the 
state dynamics during arbitrary time interval is in gen- 
eral not completely positive, but still it is linear. This 
situation is very closely related to experiments in process 
tomography, in which the state estimation of inputs is as 
necessary as the state estimation of the outputs, i.e. we 
indeed perform measurements in two time instants. And 
it could happen that already at the first time instant the 
system's and the channel's degrees of freedom are mutu- 
ally (although weakly) correlated. Fortunately, from the 
practical point of view, in this case the process tomog- 
raphy schemes are not affected, only the data processing 
should be different if one uses some advanced statistical 
tools. Wc think that this framework provides a physically 
reasonable description for the existence of noncompletely 
positive linear maps. We have shown that in a strict 
sense not all linear noncompletely positive maps can be 



7 



indeed realized within such model. The question of the 
characterization of "accessible" maps within this model 
is interesting and very important, but unfortunately we 
do not know the answer yet. 

We have argued that the problem of initial correlations 
is not a problem of quantum dynamics, but rather of 
quantum kinematics. In other words, the process tomog- 
raphy always describes a relation between the prepara- 
tors and the channel. For different sets of preparators 
physically the same channel could be described by dif- 
ferent dynamical maps. Only in very specific cases of 
preparators the channel is represented as a completely 
positive tracepreserving linear map. Fortunately, this is 
the case that usually holds in labs, or better to say, we 
are aiming to hold in our labs. We have described the 
method how to test the preparators using the calibrated 
quantum channels. Good preparator devices are crucial 



for the successful development of quantum information 
processing. The presented analysis is very far from be- 
ing complete and a deeper investigation on the system- 
environment correlations effects on quantum dynamics 
and experiments is needed. 
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